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I Introduction 



A relatively light (12-16 GeV) gluino g, along with a lighter (2-5.5 GeV) bottom 
squark b, has been proposed [1] to explain the excess of the cross section for bottom 
quark production at hadron colliders. The b squark is assumed to be a mixture of bi 
and bji, the superparterners of bi and &r. Other supersymmetric (SUSY) particles, 
except the other bottom squark b' and one of the top squarks, are assumed to be 
sufficiently heavy. The masses of b' and the light top squark t are constrained by the 
electroweak data to be below 180 GeV and 98 GeV, respectively [2]. We follow the 
convention in Ref. [1] to define 

b \ ( cos0 b sm0 b \ ( ~b R \ ( , 

V J [ -sm0~ b cos0 b )\b L ) ■ { ] 

The introduction of these new particles gives rise to new interactions in various pro- 
cesses. For example, the total decay width of the T is raised since the decay T — ► bb 
[3] is now permitted; the decay width of the Z boson is also changed [4, 5]. As a 
result, the extraction of the strong coupling constant a s at these two mass scales 
will be affected. By contributing to the /3-function, these SUSY particles slow down 
the evolution of a s with energy scale [6]. The situation has recently been studied 
in detail by Chiang et al. [7] and no clear-cut decision can be made in favor of ei- 
ther the Standard Model evolution or the evolution in the light gluino/light bottom 
squark scenario. The partial decay width F(Z — > gg) remains a key quantity to be 
determined. A better evaluation of Y{Z — > gg), among other things, can improve our 
understanding of the effect of these new particles on the electroweak measurables at 
the Z pole and hence the determination of a s [Mz) in the scenario. 

To validate the proposition of these new particles, direct searches for light gluinos 
and light bottom squarks at e~e + colliders will definitely play a key role. An analysis 
has been presented recently by Berge and Klasen [8] of gluino pair production at 
linear e~e + colliders. However, they only considered the mass range m~ g > 200 GeV. 
Production of light gluino pairs was studied by Ref. [9] and its updated version [10]. 
However, a chirally-mixed light bottom squark was not included in either of these 
calculations. The decay of on-shell Z bosons into gluino pairs was first discussed 
in Refs. [11, 12], but no chiral mixing between squarks was allowed. Djouadi and 
Drees [13] took into account chiral mixing and computed an explicit expression for 
Y(Z — > gg). However, they neglected the gluino mass and required all squark mass 
eigenstates to be heavier than Mz/2. Production of light gluinos at pp colliders was 
considered by Terekhov and Clavelli [14] but without inclusion of the light bottom 
squark either. Therefore, an analysis of light gluino production in the presence of a 
light bottom squark will be very useful for gluino searches. 

Although all the previous calculations agree with each other that the branching 
ratio of Z — > gg is less than (9(10 -3 ), they differ in some important features of 
the process. At the one-loop level the decay occurs through two types of diagrams; 
see Fig. 1. In type (a) diagrams the Z couples to a pair of quarks and a squark is 
exchanged during the process, while in type (b) diagrams the Z couples to two squarks 
and a quark is exchanged. Refs. [8, 11] disagree with Refs. [10, 13] in the relative 
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sign between the two types of diagrams. Considering only non-mixed chiral squarks 
of equal mass, Kane and Rolnick [12] claimed that the amplitude of the process is 
identically zero when m q = rriq is satisfied for each supersymmetric pair, even if weak 
isospin is broken so that, for example, m t 7^ m^. However, other references [8, 10, 13] 
state that for the contribution of quarks and squarks of a given generation to vanish, 
we must have both mass degeneracy in the quark isospin doublet (e.g., = m u ) 
and mass degeneracy in the squark isospin doublets (e.g., mg = rag, = m„ = mp). 
There are also two contradictory opinions with regard to cancellation of ultraviolet 
divergences. Refs. [8, 11, 12] asserted that ultraviolet singularities cancel separately 
for each weak isospin partner, while Djouadi and Drees [13] found that the amplitude 
is finite only after summing over a complete isodoublet. This discrepancy is essentially 
related to the relative sign between diagrams (a) and (b) in Fig. I. The divergent 
parts of the two diagrams must have opposite signs for them to cancel separately for 
each isospin partner. There may not be any constraint on the relative sign between 
(a) and (b) for divergences to cancel within an isodoublet, since the two members 
in an isodoublet have opposite I3 (the third component of the weak isospin) and the 
divergences are generally proportional to I 3 . 

A full calculation of Y{Z — > gg) involves evaluation of the Feynman diagrams in 
Fig. 1, with the cut (s) quark lines connected. To get a meaningful result, one has to 
deal with difficult one-loop integrals and remove singularities due to on-shell particles. 
In this paper we try to provide a different approach to solving the above-mentioned 
discrepancies. Since 2mb < Mz and 2m^ < Mz in the proposed scenario, the decay 
amplitude has an imaginary part which is finite and can be calculated in an easier 
way. It is likely that the imaginary part can provide some useful information on the 
full amplitude. Similar situations arise in the Kg-K^ mass difference and the decay 
Kl fi + fi~ [15]. In each case the high-momentum components of the loop diagrams 
are suppressed (through the presence of the charmed quark [16]), leaving the low- 
mass on-shell states (nn or 77, respectively) to provide a good estimate of the matrix 
element. In the light gluino and light bottom squark scenario, the decay width for 
Z — > gg usually turns out to be only a few times larger than the contribution from 
the imaginary part alone. 

This paper is organized as follows. Section II establishes the unitarity relation of 
the M. matrix elements. Amplitudes of the cut diagrams are calculated in Section III 
and the results are listed in the Appendix. The lower bound based on the imaginary 
part of the decay amplitude for Z — > gg is presented in Section IV. Implications of 
the imaginary part for the full amplitude are discussed in Sections V. Implications for 
gluino searches and running of a s are discussed in Sections VI and VII, respectively. 
Section VIII summarizes. 

II Unitarity relation 

Let us first review the decays Kl,s l l + considered in Ref. [15]. As is the case 
with Z — > gg, both decays are forbidden at the tree level. However, they can occur 
through a two-photon (77) intermediate state. Other intermediate states such as nnj 
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Z(P) f 



b(p 2 ) 
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g(k 2 ) 



Figure 1: Cut Feynman diagrams for Z — > gg: (a) Z — > (66)* — > gg, (b) Z — > (66)* — > 
gg. Similar diagrams with g(ki) <-> g^) are not shown but should be included in 
the calculation with an overall minus sign. 



and 37r are much less important. As a consequence of the unitarity of the S'-matrix 
(S^S = (1 + zT)t(1 + zT) = 1), the T-matrix element between the initial state K^ s 
and the final state satisfies the following relation 

1 r, 



Im 



(ri + \T\K L>s ) 



(l-l + \T^T\K LtS ) 



(2) 



where Im denotes the imaginary part. If we only consider the two-photon intermediate 
state, then 

d 3 k d 3 k' 



{l-l + \T^T\K L , s ) = £/ 



1 1 



(2tt) 3 (2tt) 3 2£ IE' 

x( 7 (fc,e) 7 / (A ; / ,e / )|T|^ Li5 ) 



(3) 



where {^(k, 6)7' (A/, e')) is a real two-photon state with k and k', e and e' specifying 
the 4-momenta and 4-polarizations, respectively. Since the T-matrix elements can be 
expressed as the invariant A4 matrix elements multiplied by 4-momentum-conserving 
5-functions, Eqs. (2) and (3) combine to give 

d 3 k d 3 k' 1 1 



Im 



M(K L , s ^ri- 



E 



(2tt) 3 (2it) 3 2E2E 



^7(^,6)7'^', e')\M\K LiS ) 



( 1 (k,e) 1 '(k',e')\M\ri + r(27r) i 5^(p- h - k 2 ), (4) 

times an overall 5^\p — p\ — P2), with p, p\ and p 2 being the 4-momenta of Kl,s, l~ 
and / + , respectively. It is expected that the real part of the amplitude is roughly of 
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the same order as the imaginary part, so that the actual decay width will exceed the 
lower bound based on the imaginary part by only a small factor. 

Quite similarly, the imaginary part of the invariant matrix element M.(Z — > gg) 
at the one-loop level can be written as 

Im [M(Z -> gg)] = - £ / dU f M(Z -> f)M*(~gg - mn) 4 S^(p - J», (5) 

Z / J *=1 

where the sum runs over all possible intermediate on-shell states / and dUf = 
lir=i ^^"2^" wr kh ri/ being the numbers of particles in state / and pi being the 
3-momenta of the particles. Since b is the lightest supersymmetric particle in the 
scenario and all other supersymmetric particles (except g) are expected to be heavier 
than M z /2, we only need to consider the cases where / is bb and bb. The inte- 
gral over the phase space II j can be simplified to the integral over the solid angle 
Q. In the case where the intermediate state is two particles with equal masses, we 
have / dUf(2TT) 4 5^(p — J2i=iPi) — 3^2 / dVt, where v is the velocity of the on-shell 
intermediate particles. 



Ill Amplitudes of the cut diagrams 

We adopt the spinor convention of Peskin and Schroeder [17], in which the metric 
tensor g^ u = Diagonal(l, —1, — 1, —1) and 

✓-(;;). ^f- 1 ?) - v-(j^).«-i.« «q 

where a 1 are the Pauli matrices. The uncrossed cut Feynman diagrams that contribute 
to the imaginary part of the full amplitude are shown in Fig. 1. The crossed diagrams 
with g(ki) <-> 5(^2) are not shown but should also be included in the calculation. In 
the center-of-mass frame of the Z boson, the 4-momenta of the final gluinos are 
fci = (E,k) and k 2 = (E, -k), where E = M z /2 and k = (0,0, |k|). Suppose k 
is along the z-axis and the polarizations of the Z are quantized along this axis, i.e., 
& = (0, 1, ±i, 0)/v^2 or (0, 0, 0, 1), corresponding to helicities A = ±1 or respectively. 
The 4-momenta of the intermediate bottom quarks are p\ = (E, p) and p 2 = (E, — p), 
with p = |p| (sin 9 cos 4>, sin 9 sin 0, cos 9). The 4-momenta of the intermediate bottom 
squarks are p\ = (E,p) and p 2 = (E,—p), with p = |p| (sin ^ cos0, sin 9 sin0, cos#). 
The Feynman rules for the Majorana fields are given in a representation independent 
way in [18]. 

The M. matrix element for Z — > bb is 

M{z^bb) = --^u{p 1 )Ap){glP L + g b R P R )v{p 2 )^ , (7) 

Z COS uw 

where i = e ■ 7, g\ = g b v + g\ = | sin 2 9 W - 1, g b R = g b v - g\ = § sin 2 9 W , P L = 

Pr = ^2~, ^ is a Kronecker delta in the quark color indices and p = p\ + p 2 is the 
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4-momentum of the Z. The Dirac spinors u(p\) and v (p 2 ) can be written as 

u l (pi) = 

v l (p 2 ) = 



v ] (p 2 ) = 



E-\ P \e 

e-\pW 



'e+\ 




'e-\ 




/e+\p\v 



-JE-\p\ V 



-^/e+\pW 

where the arrows j and j denote spin up and spin down along p, respectively; 



(8) 



v 



2 

e 1 ^ sin • 

— sin I 

e** cos • 



-e ^ sin 



77 



COS 2 



sin 2 



We then have 
_M(Z^& T 6 T ) 

M(Z^6 T ^) 
M(Z^& T ) 



= (0, i sin + cos cos 6>, — i cos + sin cos 6>, — sin 9) • e(p) 



x ■ 



9w 



(E-\p\)g b L + (E+\p\)g b R 



5 13 



2 cos 9w 

(0, — i sin + cos cos i cos + sin cos — sin 9) • e(p) 



x- 



2 cos 6>vk 



(E+\p\)g b L + (E-\p\)g b R 



5 ij 



gwm b 



2 cos 9 W 



gw m b ^ 
2 cos 9w 



g b L (—\, sin 6 1 cos 0, sin 9 sin 0, cos( 



+g^(l,sin#cos0, sin sin 0, cos 0) • e(p)5 tj 
g b L (l, sin 9 cos 0, sin sin 0, cos 0) 

+<7fl(— 1, sin 9 cos 0, sin 9 sin 0, cos#) • e(p)S l - 



For sin 2 6^ = 0.2311 and without top quark corrections, the partial decay width for 
Z to decay into massless bb is then ^^[(g b L ) 2 + (g b R ) 2 } = 368 MeV. 

Now we consider bb — > gg via exchange of a b or 6', the matrix element for which 
is denoted Ai(bb — > gg) or Ai'(bb — > gg), respectively. We have Ai(bb — > gg) = 
_> ^) + A^( 2 )(65 -> gg), with ' 



^(D (6 ^gg) = -2g 2 



M^&^gg) = -2g 2 



31 



(Pi - k 



1) 
(t a t b ] 



m 



(Pi - k 2 ) 2 - m 



>u(ki) (P L sin 9~ b - P R cos 9~ b ) u(pi) 
u(P2) (PflSin^ - P L cos9~ b )d(k 2 ), (9) 
^(A^CT 1 (P £ sin 9~ h - P R cos 9~ b ) u( Pl ) 



v{p 2 ) (P R sm9 b - P L cos%) Cu'ih), (10) 



where the superscript (1) denotes the uncrossed diagram and (2) the crossed diagram; 
a, b and i,j are the color indices of the gluinos and the quarks, respectively; t a are 
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the fundamental representation matrices of SU (3); C = ry°7 2 is the charge conjugate 
matrix. It can be easily verified that u(p,s) = Cv T (p,s) and v(p,s) = Cu T (p,s), 
where T means "transpose". The Majorana spinors u(ki) and v(k 2 ) also satisfy these 
relations [18]. Thus we can immediately write 




^)^^ +|k|C - 



y/E - |k|C- 



-^E+\k\(- 

-y /E-\k \C + 
E+\k\( + 



u l (k 2 ) 
v^(ki) 



v l (h) = I V ^i^t 1 , (ID 



-jE-\k\Q 



with ( + = I n ] and C- = ( 1 ) • Here the arrows j and | denote spin up and 

spin down along k (i.e., the z-axis), respectively. Since v T (k 2 )C~ 1 = —u{k 2 ) and 
Cu (ki) = v(ki), Eq. (10) can alternatively be obtained from Eq. (9) by interchanging 
ki and k 2 and adding an overall minus sign. The helicities of the final gluinos are 
determined by A, the initial helicity of the Z. For A = 1, both gluinos have spin up 
in the z-direction, while for A = — 1, both have spin down in the z-direction. For 
A = 0, one of them has spin up and the other has spin down in the z-direction. One 
expects \boa.JA(Z^ — > g^g^)\ = \lmA4(Z^ — > g^g^)\, because the two processes are 
related by mirror symmetry. One also expects lmM(Z^ — > g^g^) = 0, because these 
final gluinos have the same helicities and should therefore be excluded by the Pauli 
principle. The matrix element Ai'(bb — > gg) can be obtained from Ai(bb — > gg) by 
replacing mjj, sin^ and cos^ with my, cos^ and — sin^, respectively. 

Now consider the diagram in Fig. 1 (b) and a similar diagram with g(k\) <-> g(k 2 ), 
where the intermediate state is a pair of scalar quarks (b and b). The tree-level 
Zbb coupling is proportional to g\ sin 2 9~ b + g b R cos 2 0^, so a mixing angle of Q~ h = 

arcsin \J2 sin 2 0~w/3 ~ 23° or 157° will make it vanish. A weak Zbb coupling is assumed 
[1] to satisfy the tight constraints imposed by precision measurements at the Z peak. 
Consequently the contribution of the bb intermediate state to lmM(Z — >■ gg) should 
also be small. However, to see how the two types of diagrams shown in Fig. 1 interfere 
with each other, we take 9~ b to be a free parameter. For the first part of the cut diagram 
[Fig. 1 (b)], we have 

M (Z - bb) = [g b L sm 2 9~ b + g b R cos 2 9 b ] (p, - p 2 )%(p)^, (12) 

A COS t>w L J 

where % and j are the squark color indices. The sign discrepancy mentioned in the 
Introduction can be traced to the relative sign between Eqs. (7) and (12) [8]. The 
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current sign in Eq. (12) is consistent with the Feynman rules in Ref. [18]. We will 
argue in favor of this sign from another point of view in Section V. For the other part 
of the diagram, 

M^fb^gg) = 2g 2 ^(A^CT 1 [P L sin0 6 - P fl cos0 6 ] 

(Pi — mJ — m b 

ft l +m b )[P R smd b -P L cosd b }Cl T {k l ) (13) 

M {2 \bb^gg) = 2g 2 ( *°* ^ ^(h) [P L sin 6~ b - P R cos B b ] 

_ ^ + m& ) sin 9- b - P L cos g ] v(k 2 ), (14) 

where (1) denotes the uncrossed diagram and (2) the crossed diagram. The relevant 
matrix elements for A = 1 are presented in the Appendix. 



IV Lower bound on Y{Z 

Now we are ready to put things together and obtain a lower bound on T(Z — > gg). 
First we consider an extreme case with m b = m b = m g = and my = oo. In 
this limit, the product M(Z — > f)M(f — > gg) has an angular dependence of either 
(l + cos#) or (1 — cos0). However, the cos 9 term does not contribute to the imaginary 
part of the full amplitude, because integrating it over the solid angle Q gives zero. 
Note that tr(t a t b ) = tr(t b t a ) = 5 ab /2. The only nonvanishing amplitudes are 

M(Z^bb)*M(bb^~g^f) = 5 ab r w (g b A - g\ cos20 g ) 

M{Z^bb)*M{bb^~g^) = -S ab r w (g b A -g b v cos29 b ) 

M(Z^bb)*M(bb^g 1 g 1 ) = 5 ab r w (g b v - g b A cos29 b )cos29 b 

M(Z l ^bb)*M(bb^g l g l ) = -5 ab r w (g b v - g b A cos 29 b ) cos 29 b , 

where the asterisks imply that we have integrated over the phase spaces and summed 
over all intermediate helicity states, r w = J^^sOw " From ^ ne aDove equations we can 
see that the imaginary parts of diagrams (a) and (b) in Fig. 1 interfere destructively 
if b is more left-handed (45° < 9 b < 135°) or dominantly right-handed (9 b < 23° or 
6~ b > 157°); the contribution of diagram (b) remains negligible in the neighborhood 
of the decoupling angle (23° or 157°). The imaginary parts of the amplitudes are 

lmM(Z^ -> # T # T ) = -lmM(Z l -> g l g l ) (15) 
= 5 a V w (^-^)sin 2 ^cos 2 ^/(87r) . 

The relation (15) also holds when all the particles have a finite mass. The lower 
bound on YiZ — > gg) in the limit mj = m b = m~ g = and my = oo can be expressed 
as a ratio 

T(Z^gg) l[lmM(Z^g^g^)] 2 + [lmM(Z^g^)} 2 



T(Z ->bb) ~ 2 \hnM(Zl -> MM)] 2 + [ImM(Zi -> MM)]< 
a 2 a {g b L -g b R ) 2 ^e- b co^e- b 
6 (g b L ) 2 + (g b R ) 2 



(16) 
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The factor of 1/2 comes in because the final gluinos are identical. Taking Y(Z — > 
bb) = 368 MeV, we plot the lower bound on the decay width Y(Z — > gg) as a function 
of the bottom squark mixing angle 9 b in Fig. 2 (dotted curve). When all the masses 
are finite, we can no longer ignore the cos 9 terms, because the denominators of the 
propagators are no longer of the form ~ (1 ± cos#), which previously cancelled with 
the same factors in the numerators of the amplitudes and gave only linear terms in 
cos#. However, it is still not hard to perform the integration over the angles. Define 

t ( \ 1 r (i ± cose) 2 . 

/±(W) = 2 Jo x 2 + y 2 T# + 2xy cos 9 0d ° 

rn sin 2 9 

h{x ^ z) = Jo X i + y 2 + Z i + 2 X ycos9 Sined6 ' 
and let c± = I±(v b ,v g ,r~ b ), c' ± = I±(v b ,v g ,ry), c = I (v b ,v g ,r b ), c' = I (v b ,v g ,ry), 



c = I (v~ b ,v g ,r b ), where r { = 2mi/M z (i = b,b,b',g), v { = yjl - r 2 is the "velocity'^ 
of an on-shell particle i (i — b, b, g) . The lower bound can then be written as 
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nz - gg) > ^^(^ + + A ^)\ » ( 18 ) 

where, up to a common factor of proportionality, A±v b and A[v b are the imaginary 
parts of the amplitudes for Z — > (bb)* — > gg via exchange of a b and b', respectively; 
A 2 v b is the imaginary part of the amplitude for Z — > (bb)* — * gg. We have 

A 1 = c 1 (g b A -g b v cos29 b )-g b A (c 2 + c 3 sm29 b ) 
A = c[{g b A + g b v cos29 b )-g b A (c' 2 -c' 3 sm29~ b ) 
A 2 = c v~v g (gy - g\ cos 29 b ) cos 29 b , 

where c x = c_(v g - v b ) + c + (v g + v b ) + c r^, c 2 = (c_ + c + + c )r 2 v g , c 3 = 
(c+ — C-)v b r b r g ] c[, c' 2 and c' 3 are defined similarly, with c± and Co all primed. All 
these quantities only depend on the masses. As m b , m b and m g go to zero and my goes 
to infinity, C\ — > 1, c v?v g — > 1, c 2 — > 0, c 3 — > 0, c' x — > 0, c 2 — > and c 3 — > 0. Eq. (16) 
is thus recovered. As far as the imaginary part of the amplitude is concerned, Eq. (18) 
agrees with Ref. [13] in the limit m g = except for the mentioned sign discrepancy. 
Compared to Ref. [13], this unitarity calculation not only takes into account a nonzero 
gluino mass but also is much simpler. Firstly, there are no singularities to remove. 
Secondly, one only has to evaluate a few elementary integrals in Eq. (17) instead of 
the much more difficult one-loop two- and three-point functions. 

The lower bound is plotted in Fig. 2 (solid curve) as a function of 9 b for a specific 
set of values for the masses: m b = 4.1 GeV, m b = 4.5 GeV, my = 170 GeV, m g = 15 
GeV. The small peak around 90° disappears if the gluino has a small mass (e.g., 
m g < 7 GeV) such that the sign of the sum A\v b + A\v b + A 2 v b does not change over 
0° < 9 b < 180°. If g were massless, the lower bound as a function of 0° < 9 b < 180° 
would have two identical peaks (see the dashed curve in Fig. 3). Due to destructive 
interference between diagrams (a) and (b) for 45° < 9 b < 135°, 9 b < 23° and 9 b > 157°, 
the lower bound is smaller than the contribution from the bb intermediate state alone 
(Fig. 2, dashed curve) in these ranges. 
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Figure 2: Lower bound on F(Z — > gg) as a function of the bottom squark mixing 
angle 6~ b . Solid curve: m b — 4.1 GeV, m~ h = 4.5 GeV, m- g = 15 GeV, my = 170 GeV; 
dashed curve: only the contribution from the bb intermediate state is included with 
the same set of masses; dotted curve: m b = m~ h = m-g = and my = oo. 
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V Implications for the full decay width F(Z 



Let us first examine under what conditions the decay amplitude vanishes. We will 
see that Kane and Rolnick's claim [12] contradicts our unitarity calculation. Under 
their assumption, m q = m^ = m q and 9q = 0, so that q = c' { (i = 1,2,3) and 
Co = c' = Io(vi>,Vg,r b ) = ±c . We regard the relative sign (±) between diagrams (a) 
and (b) as undecided and only consider the contribution of the third generation, i.e., 
q = t,b. Define A' 2 = — c' VyVg(g v + g A cos 29~ b ) cos 29~ b . Then the imaginary part of the 
amplitude is proportional to A\ + A\+A2 + A 2 = 2[(c + + C-^fCo)v b Vg + (c + — c^)vb]g b A , 
which is nonzero for either sign. Therefore, the amplitude does not vanish under Kane 
and Rolnick's conditions. For the imaginary part to be zero, we must have m t = m b 
and mi = m t -, = m~ h = m b , and similar mass degeneracies in the other two generations. 
(Of course, the imaginary part of the amplitude automatically vanishes if all quark 
and squark mass eigenstates are heavier than Mz/2.) This is true whether the relative 
sign between diagrams (a) and (b) is + or — and is consistent with Refs. [8, 10, 13]. 

Let us now investigate whether the imaginary part of the amplitude computed in 
the previous section can give us some hint how loop divergences cancel and reasonably 
small decay widths for Z — > gg can be obtained. Without actually calculating the full 
amplitude, we should be able to recover part of it from the imaginary part. For sim- 
plicity we take m~ g = and consider only diagram (b) in Fig. 1. The imaginary part of 
the amplitude for this diagram is proportional to A2V1 = 5ov~(gy — g\ cos 29 b ) cos 29 b . 
We have 



7T _ o 

2 CoV ' b 



1 t?t| — ml 






v~ b ir + 


2 s 



2 , K-™') 2 " 

m b H 



n h (ml-m 2 b )/s-(l-v~ b )/2 
s (ml-m 2 b )/s-(l + v b )/2 : 



where s = Mf . For m b < M z /2, one can check that 
lmB (s,m b ,m- b ) 



V b 7T 



T r( \ 1 (ml-ml)/s-(l-v b )/2 

lmC (s,m b ,m- b ,m b ) = -log _ _ (i - — 



where B and Co are the scalar one-loop two- and three-point functions [13, 19, 20, 21], 
respectively. If we define 



A(m 1 ,m 2 ,m 3 ) 



1 (mf — mf) + [rn\ — m 
2 



+ 



2s 

2 (m\-ml)(ml-ml) 
rn 3 -1- 



B (s,m 1 ,m 2 ) 
C (s,mi,m 2 ,m 3 ) 



then by analyticity and symmetry, the following terms should be part of the full 
amplitude contributed by type (b) diagrams in Fig. 1, 



A ( m h m h m b)(g b v - g A cos20 b) cos2d b - A ( m b>i m b>i m b)(g b v + g b A cos20 b) cos 2e i 

-A(m b , m b ,,m b )2g b A sin 2 29~ b + A{m h m h m t ) (g v - g\ cos 20f) cos 29~ t 

-A(77ip, ?7ip, m t )(g v + g\ cos 2%) cos 29~ t - A(m h m- t ,, m t )2g\ sin 2 29~ t . (19) 
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Ml = 





m 



Here the sin 2 26 q terms can be obtained by repeating the unitarity calculation with 

/ = bb'. Alternatively, they can easily be guessed if we note that the above terms 
should sum up to zero for m t = m b and m^ = my = m~ b = my. Those terms in Eq. 
(19) are expected to be the only ones relevant to our argument [13]. 

The top and bottom squark masses and mixing angles are determined by the 
following mass matrices [13] 

™ 2 i L + ™?t + g\s cos 2/3/2 -m t (A t + (j, cot (3) \ 

-m t (A t + /i cot (3) m 2 iR + m\ - g l R s cos 2/3/2 J ' ' 

$ hL +m 2 b + ^scos 2(3/2 -m b (A b + //tan/3) \ 
-m b (A b + fj, tan (3) m\ R + m\ - cos 2/3/2 ) ' ^ > 

where m 2 ^m 2 ^m 2 h ^m 2 h ^ are soft SUSY breaking masses; tan/3 is the ratio of the 
vacuum expectation values of the two neutral Higgs fields in the MSSM; A t) A b de- 
note the trilinear Higgs-stop, -sbottom couplings, respectively; and /t is the Higgsino 
mass parameter. 577(2) gauge invariance leads to m~^ = m~ . It is reasonable to 
assume that we should get a sensible decay width Y(Z — > gg) for each specific set of 
parameters. In particular, we can choose very large values for m^, mp , m~ b , my and 
should find a tiny T(Z — > gg) for some mixing angles Qi and 9~ b . For simplicity we 
assume A t + ji cot (3 = A b + ji tan (3 = so that there is no L-i? squark mixing, and 

m\ h + g\s cos 2/3 /2 = m\ R - g R s cos 2/3 /2 
m 2 L + cos 2/3/2 = m~ fl - cos 2/3/2 

so that m 4 ~ = , m^, = m^,. The mass difference between the top and bottom squarks 
is then 

m t ~ m l — J 171 } + m t + 9l s cos 2/3/2 — ^m 2 + m 2 + g£s cos 2/3/2 ps , 



m f - 



with <5m 2 = [2(m 2 — m 2 ) + (g^ — g^)s cos 2/3]/4. The sum in Eq. (19) becomes simply 

A ( m b i m h m b) ~ A ( m h m b m t) = SA . (22) 
In the heavy squark limit, 



2 

TfL ~ 

B (s,m q ,m q ) = A + log— I 



2 2 

m^ m q 

C {s,m q ,m q ,m q ) = + — ^log-f 

m| — m q (m| — m q ) A m q 



where A denotes the ultraviolet divergent part of B Q and v is the renormalization 
scale. We can see that the B and C terms in A(m q ,m q ,m q ) vary as — (to- — 
m 2 ) [log m~/z/ 2 ]/s and — (m~ — m 2 )/s, respectively. So the leading term in SA comes 
from B and varies as 

[(mf - to 2 ) logm? - (to 2 - m 2 b ) log mf]/s ~ [(^ - g b L ) cos 2/3 log m 2 J /2 . 
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Thus the sum in Eq. (19) is logarithmically divergent as m~ tL goes to infinity. There- 
fore, we cannot get reasonable decay widths if this divergence is supposed to cancel 
within an isodoublet. So we do not agree with Djouadi and Drees' claim [13] that 
one can get meaningful results by summing over a complete isodoublet. The only 
other way out is for the divergence to cancel for each weak isospin partner. Note 
that the ultraviolet divergent parts of the B terms in Eq. (19) not only can cancel 
within an isodoublet, but also can cancel for each weak isospin partner. Indeed, if the 
relative sign between diagrams (a) and (b) in Ref. [13] is reversed, one finds that both 
ultraviolet divergences and the divergences in the heavy squark limit (as discussed 
above) cancel separately for the top and bottom sectors. To understand how the 
latter divergences cancel, repeat the same reasoning as that leading to Eq. (19) for 
type (a) diagrams in Fig. 1 and observe that those diagrams should contribute the 
following terms, 

± \C(m b ,m b ,m b )(g b A - g b v cos26~ b ) + C(m b ,m b ,m b ,)(g A + g b v cos 2%) 

C(m t , m t , mi) (g A - g l v cos 20 1) + C(m t ,m t ,mi,) (g l A + g l v cos 20 t ~)] , (23) 

where again we keep an undecided relative sign and 



C(m 1 ,m l ,m 2 ) 



9 (ml - ml) 2 
m\ + — — 



C (s,m 1 ,m 1 ,m 2 ) 



For m t - = m t -, and m~ b = my, Eq. (23) becomes 

±[C(m t , m t , mi) - C(m b , m b , m b )] = ±5C . (24) 
One can check that when m q becomes large, 

S 771- 

ReC (s, m q ,m q , mi) ~ ~ 9 log — . 

m~ — m z q s 

Thus the leading term in Re<5C is exactly the same as the leading term in 5 A. The 
second-to-leading and higher order terms turn out to be finite after summing over 
q — t,b (for either relative sign) or over diagrams (a) and (b) (for only one of the 
signs). If we have chosen a correct relative sign, the divergences in type (b) diagrams 
should cancel exactly with those in type (a) diagrams separately for q = t,b. Otherwise 
the divergences will add up and the total amplitude will be logarithmically divergent 
as the squark masses go to infinity. We find that the sign in the current note rather 
than that in Ref. [13] is favored. 

Our numerical analysis not only verifies the above argument but also shows that it 
works even if arbitrary L-R squark mixing is allowed and the above constraints on the 
parameters are relaxed. Using the same formula in Ref. [13] but with the sign flipped 
3 , we find that for m- g ~ 0, m b < (9(30) GeV, my = 0(150) GeV, m~ t = O(90) GeV, 
mil = (9(300) GeV, T(Z — > gg) is typically of order 0.1 MeV if b is not dominantly 

3 The formula fails to produce sensible decay widths in the heavy squark limit if the sign is not 
reversed, as argued in the text. 
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left-handed and only of order 0.01 MeV or less if b is dominantly left-handed, i.e., 
6~ b pa 90°; and that the top squark mixing angle Of has little effect on the decay width. 
If all squark masses are greater than M z /2, we confirm the statement in Ref. [13] 
that YiZ — > gg) depends weakly on the details of L-R squark mixing and find that 
it is of order 0.1 MeV for a wide range of MSSM parameters. Thus inclusion of a 
dominantly left-handed light b would reduce the decay width YiZ — > gg), as is the 
case for the unitarity lower bound plotted in Fig. 2. A relatively heavy or not very 
left-handed b will not change the decay width by much. In Fig. 3 we plot the full 
decay widths as well as the corresponding unitarity lower bounds for m g = and two 
sets of squark masses. When all the squark mass eigenstates are heavier than Mz/2, 
the full width is about an order of magnitude larger than the lower bound. When 
one of the bottom squarks is light (< O(30) GeV), the shape of the full width as 
a function of 9~ b is similar to the lower bound and generally only a few times higher 
than the latter. Although an expression for YiZ — > gg) is not available when trig ^ 0, 
we expect the shape of YiZ — > gg) to be also similar to and only a few times higher 
than the unitarity lower bound plotted in Fig. 2 (solid curve) if the gluino has a 
mass around 15 GeV and the light bottom squark has a mass around 4.5 GeV. So 
T{Z — > gg) should be of order 0.1 MeV in the light gluino and light bottom squark 
scenario. 



VI Implications for gluino searches in Z decays 

Aside from Z — > gg, there exist three other gluino-producing Z decays, Z — > bbg, 
Z — > bbg and Z — > qqgg. The first two processes are ~ aa s at the tree level and have 
a combined decay width of 1.9 - 5.9 MeV depending on the sign of sin 26~ b [22]. The 
third process is ~ aa 2 s and its decay width is calculated in a model-independent way 
to be 0.75 - 0.21 MeV for m g =12-16 GeV [23]. A recent analysis [24] shows that 
Y(Z — > bbgg) can be enhanced by 10% - 60% due to additional "sbottom splitting" 
diagrams. This will raise YiZ — > qqgg) by 0(0.01) MeV. The new SUSY particles 
do not always contribute positively to the Z width, however. Cao et al. [5] and S.w. 
Baek [25] showed that the decay width YiZ — > bb) can be reduced by 2 - 8 MeV. By 
fine-tuning the parameters in the light gluino and light bottom squark scenario, all 
the electroweak measurables (Yz, Yh a d(Z), Rb, R c ) at the Z pole can be still within 
the la - bounds of the experimental values. Thus, existence of the new particles can 
only be verified through direct searches for gluinos or bottom squarks. The light 
bottom squark is assumed to be long-lived at the collider scale or to decay promptly 
to light hadrons in this scenario. In either case, it forms a hadronic jet within the 
detector due to its color charge, g decays exclusively to bb or bb and becomes two 
hadronic jets. The smallness of YiZ — > gg) implies the insignificance of Z — > gg in 
gluino searches. Searches for signals of Z — > bbg + bbg and Z — > qqgg will be expected 
to play a pivotal role. 
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m =4.5 GeV 

b 




Bottom Squark Mixing Angle (degrees) 



Figure 3: T(Z — > gg) as a function of the bottom squark mixing angle 9~ b . We take 
6~ t = 50°, m b = 4.1 GeV, my = 170 GeV, m t = 17 A GeV, m~ t = 95 GeV, m~ t , = 300 
GeV, mj = 0. Solid (full width) and dashed (unitarity lower bound) curves: = 4.5 
GeV; dash-dotted (full width) and dotted curves (unitarity lower bound): = 100 
GeV. 
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VII Implications for running of a s 



Both the light gluino g and the light bottom squark b can change the /^-function that 
governs the energy-scale dependence ("running") of the strong coupling constant a s . 
At the two-loop level, a s (Mz) can be raised by 0.014 ± 0.001 [7] with respect to its 
standard model value if extrapolated from the mass scale m b . A natural question 
arises: are values of a s (M z ) determined from measurements at different energy scales 
still in accordance in the presence of g and 6? To answer this question, the effects of 
the new SUSY particles on measurements at different scales must be analyzed. For 
example, the hadronic width of the Z is changed in two ways: 1) the interference 
of the standard model diagrams and the diagrams with the SUSY particles in loops 
will reduce the partial width of Z — > bb; 2) the existence of the new decay channels 
Z — > bb, Z — > gg, Z — > bbg/bbg and Z — > qqgg will raise the hadronic width. The 
bottom squark mixing angle 9~ b is chosen to be near 23° or 157° so that Y(Z — > bb) 
is suppressed 4 . Y(Z — > gg) is only of order 0.1 MeV at either of these angles. Thus 
these two channels combined will change the predicted hadronic width of the Z by 
a very small amount compared to the decrease in YiZ — > bb) and the increase in 
rhad(-/j) due to Z — > bbg/bbg and Z — > qqgg. A better determination of Y(Z — > bb), 
Y(Z — > bbg/bbg) and Y(Z — > qqgg), or a more precise measurement of Rj, (which will 
constrain the value of Y(Z — * bb) more tightly), is needed for a clear-cut decision 
in favor of either the Standard Model or the light gluino and light bottom squark 
scenario. 

VIII Summary 

We have calculated the imaginary part of the decay amplitude for Z — > gg and 
used it to analyze the full amplitude and solve some discrepancies in the literature. 
We have confirmed the argument that the decay width vanishes if both quarks and 
squarks of a given generation are degenerate in mass (the quarks and squarks in that 
generation do not need to have equal mass). We have found that both divergences 
in the heavy squark limit and ultraviolet divergences cancel for each weak isospin 
partner, as previously claimed by Refs. [11, 12]. We also favor the relative sign of 
Refs. [8, 11] between diagrams (a) and (b) in Fig. 1. Borrowing the formula for 
Y{Z — > gg) from Ref. [13] but with their relative sign between diagrams (a) and (b) 
flipped to be consistent with our calculation, we find that the decay width is of order 
0.1 MeV in the proposed light gluino and light bottom squark scenario. Compared 
with other decay processes like Z — > bbg/bbg and Z — > qqgg, Z — > gg will only play a 
moderate role in searches for gluinos and analysis of effects of the SUSY scenario on 
a s (M z ). 

4 T(Z -> bb) will be greater than 15 MeV if 40° < 6~ h < 140°. 
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Appendix: relevant M-matrix elements 

We define 

Ai i) = _ 2g 2 (t b t a h 

= -2g] 
i« = _ 2 , 2 
i (2) = -2gl 



(Pi 


-fci) 2 - 


ra- 
ft 








(Pi 


- A; 2 ) 2 - 


ra~ 






Wii 




(Pi 


- &i) 2 - 


ra 2 








(Pi 


- A: 2 ) 2 - 


ra 2 



5 



B±± = J(E±\k\)(E±\k\ 



S- b = sm9 b 
C b = cos9 b 



M. matrix elements for Z 1 * — > 66 — > ^(jf: 

■^V^ [(e-ip|)^+(e+ip|)^' 

V 2 cos 6^ L 



M(Z^ - 


► 6 T 6 T ) 










7W(Z T - 












M(Z T - 













1 + cos 9 ■ ■ 

<5 J 

2 



= -^e^E^ - S- + C 5 )(S + -5 S - 5_ + C7 s )I±|2^ 
= A«e-^(S__5 S - B ++ C b )(B__S b - B ++ C b ) 1 -±^ 

V 2 costly ^ 



1 — cos 9 



= -A^e^{B_ + S b - B + _C b ){B_ + S b - B + -C~ b ) 1 -^ 

= ~ 9 r Wmh (9 b L + 9 b R )—5- 

= -A^{B + _S b -B_ + C^B ++ S b -B__C b ) S ^ 

= A^\B__S b -B ++ C b )(B_ + S b -B + _C b f^ 
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cos 9 



M W(b l F -> ^¥) 
M^ibW -> ^ T ) 



COS 



sin # 



= -A^e" 2 "^^ - B + -Ci)(B-S- b - B ++ C h ) 



2 

sin6> 



.M matrix elements for Z 1 * — > 66 
M(Z T ^66) = 



0w 



v^cosfl 



" e P 



sin6> 



.M {1) (66 



AW e -^|p|(iL_Sf + £++C£) sintf 
-i (2) e^|p|(5 ++1 Sf + £-_Cf) sin# 
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